Non-equilibrium electronic transport and interaction in short metallic nanobridges 
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We have observed interaction effects in the differential conductance G of short, disordered metal 
bridges in a well-controlled non-equilibrium situation, where the distribution function has a double 
Fermi step. A logarithmic scaling law is found both for the temperature and for the voltage depen- 
dence of G in all samples. The absence of magnetic field dependence and the low dimensionality of 
our samples allow us to distinguish between several possible interaction effects, proposed recently 
in nanoscopic samples. The universal scaling curve is explained quantitatively by the theory of 
electron-electron interaction in diffusive metals, adapted to the present case, where the sample size 
is smaller than the thermal diffusion length. 
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I. INTRODUCTION 

Dynamically screened electron-electron interactions in 
disordered metals are known to cause singular correc- 
tions to the electronic density offStates (DOS), as ex- 
plained by Aronov and Al'tshulerllla. The correspond- 
ing zero-bias anomalies (ZBA) have been studied by tun- 
neling spectxQSCopy on wide junctions in thermodynamic 
equilibriumEta. How is the ZBA modified when a large 
current is driven through the system by a finite bias 
voltage, and the bridge is so small that there is phase 
coherence across the whole sample? In this case, no lo- 
cal equilibrium is reached, since energy relaxation does 
not occur within the sample. Instead, the distribution 
function of electrons traversing the system exhibits two 
Fermi steps related to the two different electrochemical 
potentials oLthe two leads. This fact has been claimed 
theoreticallyQij ancUhas recently been confirmed in tun- 
neling experimentsQ (where in addition the steps were 
rounded due to interactions in long wiresjLand indirectly 
by noise measurements in a similar setupu. 

We have fabricated metallic nanobridges, which al- 
low for ±he first time to study the Aronov-Al'tshuler 
anomalyEI in the stationary non-equilibrium situation 
specified above. A two-dimensional (2D), diffusive metal 
bridge, much shorter than all inelastic scattering lengths, 
is placed in good metallic contact between two thick 
reservoirs. A high lead-to-bridge aspect ratio allows 
to maintain a finite voltage drop across the bridge de- 
spite its small resistance, thus allowing the observa- 
tion of temperature-voltage scaling behavior. The study 
of a 2D bridge is further motivated by the fact that 
this geometry allows to distinguish between the Aronov- 
Al'tshuler conductance anomaly and a possible tw| 
channel Kondo effect induced by two-level system; 
which has been put forward as ther-arigin of ZBAs ob- 
served in ultrasmall point contactsty; In 3D both the 



Aronov-Al'tshuler and the two-channel Kondo anoma- 
lies show square root power- law behavior; in 2D the 
Aronov-Al'tshuler correction is logarithmic, while the 
two-channel Kondo singularity, as a local effect, is in- 
dependent of dimension. As will be discussed, the con- 
ductance through the non-equilibrium bridge can be well 
described by a two-termnal Landatter-Biittiker formula, 
generalized for interacting systemgl^l. Applying the the- 
ory of electron-electron interaction in disordered systems 
in non-equilibrium, we find that the electronic DOS is 
strongly modified by the non-equilibrium and displays 
two logarithmic singularities, corresponding to the two 
steps in the distribution function. This yields quantita- 
tive agreement with the experimental results. 




FIG. 1. SEM picture of one of the samples, which were 
produced by the shadow-evaporation technique. The thin 
Cu82Aui8 bridge (dark) is embedded between two con- 
siderably thicker Cu leads. Dark regions are the thin 
shadow-evaporation replicas of the thick leads. A thick replica 
of the bridge, which is not connected to the electrodes, is re- 
moved from the picture for clarity. 
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II. EXPERIMENT 



The sarpples were produced by the shadow-evaporation 
techniques. Metal was deposited at two different an- 
gles under UHV conditions through a PMMA mask pat- 
terned with e-beam lithography. We deposited 10 nm of 
Cuioo-2:Aua; for the bridge at one angle (-13° off the nor- 
mal direction of the surface) and 700 nm of Cu for the 
reservoirs (-1-13°) at the opposite angle without break- 
ing the vacuum (on the 10~^ mbar scale during evapo- 
ration). The transport measurements were performed in 
a ■^He/^He dilution refrigerator in a temperature range 
between T = 100 mK and T = 2.1 K. The differential re- 
sistance was measured with a LR 700 resistance bridge, 
which was capacitively coupled to a DC circuit through 
the sample. The AC excitation was less than 5 /iV. 
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TABLE I. Characteristic parameters for the four investi- 
gated nanobridges: material composition, thickness d, resis- 
tance R, diffusion coefficient D, mean free path I, and the 
amplitude A of the zero-bias anomaly at an applied magnetic 
field of -B = and B — 8.5 T, respectively (see text). In high 
magnetic field the pure Cu sample exhibits the same ZBA as 
all other samples, but in low field shows an enhanced am- 
plitude of the effect, while the scaling behavior Eqs. (1), (2) 
remains valid. This needs further investigation. 

The parameters characterizing the four different samples 
investigated in this study are summarized in table I. Each 
of the four nanobridges is L = 80 nm long, about 80 nm 
wide and their thickness d is substantially smaller than 
the lateral extension L. It is placed in good metallic 
contact between two bulk Cu leads, which are about 70 
times thicker than the bridge and extend over a large 
area of about 1 mm^ each (see Fig.||). Hence the volt- 
age applied to the sample drops only along the bridge, 
and the arising heating power is reliably conducted away 
by the leads. This was verified in a control experiment, 
where the cooling capability of the reservoirs was varied 
using a different geometry. No change of the behavior re- 
ported below was observed. Moreover, as demonstrated 
in detail below, the electronic motion is coherent over the 
entire size of the nanobridge, which means, in particular, 
that there is no heat deposition within the bridge. The 
elastic mean free path i was estimated from the residual- 
resistance ratio and consistently from the Drude formula. 
In all samples £ is comparable to the thickness d, but 
much shorter than the lateral length L, and therefore, 
the electronic density modes in the bridge obey the rules 
of 2D diffusive motion. Below we present experimental 



conductance data of sample no. (1); the results obtained 
from the other Cuioo-xAua, and Cu films are very similar. 

Raw data of the zero-bias conductance as a function 
of temperature T are shown in Fig. ||. We observe a 
logarithmic T dependence between T = 100 mK and 
T = 2 K: 

G(0,T) = G{0,To = lK) + ^-ln(r/ro) , (1) 

where A denotes the amplitude of the effect, with A = 
0.49 e'^/h for this particular bridge. Below T = 100 mK, 
small deviations from the form Eq. (1) were observed, 
which are probably due to incomplete thermalization of 
the sample at the lowest temperatures. Such logarithmic 
behavior has often been observed in thin metallic films 
and metallic naiiO|bridges and is usually attributed to 
weak localizationEJ and/or electron-electron interactiontl. 
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FIG. 2. Zero-bias conductance of the Cu82Aui8 bridge. 
The straight line is a fit according to Eq. (|l|). 

When applying a bias voltage U, small deviations from 
Ohmic behavior, i.e. a voltage-dependent differential 
conductance G = dI{U)/dU, are observed at low U, as 
shown in Fig. |^. This zero-bias anomaly (ZBA), sym- 
metrical with respect to reversal of U , becomes continu- 
ously more pronounced as T is lowered. For all voltage 
sweeps taken for various T in the range of 0.1 K to 2.1 K, 
one observes a remarkable scaling property: If one sub- 
tracts the zero-bias conductance of each sweep and plots 
G(J7, T) - G(0, T) as a function of U/T, all conductance 
sweeps collapse onto a single curve, as can be seen in 
Fig. ^. Comparison with results of several Cuioo-^Aua, 
samples and for pure Cu samples (with a mean free path 
comparable to the thickness of the wire) yields a further 
general property (G(C/,T) - G{0,T))/A is identical for 
all samples (table I), where A is the zero-bias anomaly 
amplitude defined above. We therefore propose the fol- 
lowing scaling law: 

^ ^^^U/ksT) (2) 
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FIG. 3. Top: Raw data of the voltage dependence of the 
differential conductance for different temperatures: T=0.1 K, 
0.25 K, 0.5 K, 0.75 K, 1 K, 1.5 K, 2.1 K. Bottom: The scaling 
of the same data appears, when the zero-bias conductance is 
subtracted from each curve and the data are displayed as a 
function of eU /ksT. 

where $ is a function of eU /ksT only. Other scal- 
ing laws, in particular ^/T depjendence as expected for 
the two-channel Kondo effectliy, fit very poorly to the 
data. Fig. ^ shows the function $ extracted from con- 
ductivity measurements on one sample at seven differ- 
ent temperatures as a function of y — eU/kBT on a 
semi-logarithmic scale (data of positive and negative bias 
are included). Within some scatter due to experimen- 
tal noise, all sweeps collapse onto one single curve. At 
|y| < 1, the conductance is essentially constant, which 
corresponds to thermal smearing. At \y\ > 1, the be- 
havior is roughly logarithmic, and for still larger values 
of 1 2/1, the curves deviate from scaling. The low-T data 
follow the logarithmic behavior up to y ~ 100, while 
the high-T data deviate at lower y. These deviations at 
higher energies are attributed to heating effects and/or 
additional high-energy processes like phonon excitations 



and are well known from similar experimentaiS. 

When a perpendicular magnetic field of up to -B = 8.5T 
is applied, we do not observe any systematic change of 
the conductance behavior apart from universal conduc- 
tance fluctuations (UCF). In particular, the scaling be- 
havior persists, with the amplitude A and the scaling 
function $ remaining unchanged. We expect a weak lo- 
calization conductance correction of SGwl{U — 0,T — 
0, -B = 0) w 3.2e^/?i, which is not clearly distinguishable 
from the UCFs in our small samples. However, the fact 
that the amplitude A of the ZBA as a function of applied 
voltage U (i.e. the slope of the logaritmic scaling curve) 
remains unchanged in a magnetic field of _B = 8.5T indi- 
cates that weak localization only gives a constant, voltage 
independent correction and that the logarithmic increase 
of the conductance with applied voltage is not due a to a 
loss of coherence at finite voltage and a subsequent sup- 
pression of weak localization, since otherwise the effect 
would disappear in a magnetic field. This, in turn, pro- 
vides clear evidence that the samples are coherent even 
at finite bias voltage U as expected from their small size. 
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FIG. 4. Logarithmic scaling plot according to Eq. (^ for 
the data presented in Fig. 3. Lines represent experimental 
data for various fixed T. Circles represent the theoretical 
scaling curve computed from Eqs. (^, (^o|), (^. 

Logarithmic behavior may also be caused by mag- 
netic impurities or by non-magnetic two-channel Kondo 
defectsEI above their respective Kondo temperatures Tk- 
Since an applied field of i? = 8.5 T does not modify the 
ZBA, any magnetic impurities present in the sample must 
have Tk ^ 8.5 K. However, the logarithmic behavior of 
the zero-bias conductance observed down to the lowest 
T (Fig. 1^) puts an upper bound to the Kondo temper- 
ature, Tk < 0.1 K, thus ruling out magnetic impurities 
as the origin of the ZBA. In the two-channel Kondo sce- 
nario, from poini^rpcitact spectroscopy on Cu one expects 
Tft: « 5 to 10 KEIM Hence, it is unlikely that the ZBA 
is due to two-channel Kondo defects for the same reason 
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as in the magnetic case. The assumption that there is no 
sizable number of two-channel Kondo defects present in 
our Cu82Aui8 samples is consistent with the fact that in 
previous Cu point contacts the two-channel Kondo sig- 
nal completely disappeared upon doping with 1% Au or 
moreEEl. 



m, n of the bridge. 

Using the Dyson equation for the Keldysh matrix Green 
function Qnaii-^) in terms of the exact eigenstates of the 
bridge, it is straight-forward to show that at finite bias 
and for negligible energy relaxation the distribution func- 
tion f{uj) of quasiparticles is a linear superposition of the 
distribution hmctions in the leads, 



III. THEORY 

In the following we show that the observed scaling be- 
havior of the differential conductance can be explained 
quantitatively by electron-electron interaction within a 
model calculation for the non-equilibrium DC transport 
through our small, disordered metal bridges. Since the 
diffusion time of an electron through the bridge, = 
L^/D « 1.8 ps, is short compared to the dephasing time 
and the energy relaxation time r,;, the electrons oc- 
cupy the exact eigenstates of the disordered bridge while 
traversing the system. Therefore, the DC transport is co- 
herent (zero-dimensional), even when a finite bias volt- 
age is applied as discussed above, i.e. the interacting 
eigenstates of the bridge serve as transmission channels. 
In such a situation a generalized Landauer-Biittiker ap- 
proach can be applied, where the conductance is ex- 
pressed in terms of the interacting, non-equilibrium DOS 
J\f{u!\-oi the bridge in the presence of coupling to the 
leadstll. Diffusive density modes, which exist at time 
scales down to the elastic scattering time r = i/vp ~ 
4.1 fs ^ td couple to the electronic DOS via the dy- 
namically screened Coulomb interaction and thus give a 
singular correction to the conductance as shown below. 

Because of the_good metallic contact, charging effects 
at the interfacestJ are negligible in our devices. Since 
the transition between the reservoirs and the bridge oc- 
curs abruptly on a scale small compared to the extension 
of the bridge wave functions and energy relaxation in 
the bridge is negligible, one may represent the system in 
terms of the exact eigenstates jncr) of the bridge, sepa- 
rated from the leads, and include the bridge-lead, cou- 
pling to infinite order in perturbation theoryllj^tj (as- 
sumed to be convergent). Denoting the corresponding 
creation operators for eigenstates with spin a in the left 
(right) leads and the bridge by cjjg-j 
tively, the hamiltonian of the device is H — Hl 



respec- 
are 



Hb+Hv- Here, H^n-^ = Efe efeaL(ii) fca»L(fl) fc^ 
the left (right) lead hamiltonians with spectral density 
Aka{i^) and a flat bulk DOS No = J2ka ^ka{uj)- Hy = 

'EknA^nkcLaL ka + V^kcLaR + dcscribcs mixing 
between bridge and leads via transition matrix elements 
V^f} , and Hb is the hamiltonian of the bridge, includ- 
ing random impurity scattering and Coulomb interaction. 
For later use we also introduce the effective couplings 

r^(«), and 



rL{R) * _ 



r = r^r^/(r^ + r^), which win be assumed to be as- 
sumed to be independent of energy and of the channels 



fiE) = 



1 



[ri/(")(S) + r^f(°){E + eU)] , (3) 



where f°{E) = l/(e^/'=«^ -t- 1) is the Fermi function. 
The distribution function f{E), Eq. (||), represents the 
probablility for eigenstates of the entire bridge with en- 
ergy E being occupied and, therefore, is position inde- 
pendent. It should be emphasized that this is to be 
distinguished from the local distribution function fx{E), 
which is defined as the probability of finding an elec- 
tron in a state with energy E within a small volume 
centered at position x wherein the externally applied 
fields are slowly varyingO. When, as in the present 
case, the interaction correction to the local DOS is small 
{5N^{E)/Nx{E) < 1), in a di^ered system f^{E) 
obeys a diffusive kinetic cquatior 



DV^UE)^C{{f^},x,E). 



(4) 



For negligible energy relaxation the collision integral C 
vanishes, and fx{E) is a linear superposition of Fermi 
functions f'^^^E), f'^^^E + eU/h) with coefficients de- 
pending linearly on the position x along the bridge. It 
follows directly from the definitions of f{E) and fx{E) 
via the total and the local density matrices, respectively, 
that 



fiE) 



dx 



(5) 



The local distribution fx{E) has been measured in tun- 
neling experiments of Ref. [7] , while for the description of 
the present experiment only the integrated distribution 
f{E), Eq. (|), is needed. 

An exact, Landauer-Biittiker- like expression for the 
current I{U) follows frorip_the equation of motion of the 
reservoir electron densityllil. 



dE 



r{E) ~ riE + eU) rtr [lu,g-{E% (6) 



where Q"^ is the exact retarded single-particle Green func- 
tion for the bridge in the presence of the leads, and the 
trace extends over the complete basis of bridge eigen- 
states, including spin. The problem of the DC current at 
arbitrary bias voltage is thus reduced to calculating the 
non-equilibrium DOS, M{E) = -^tr Im G^'iE). Note 
that the diffusion constant D does not appear explicitly 
in Eq. (0) , since the DC transport is zero-dimensional. 



4 



p CO 



p-q co-Q 



P CO 



FIG. 5. Leading DOS correction. Solid lines: conduction 
electron Green functions; wavy line with solid triangles: dy- 
namically screened Coulomb interaction Vq. 

The leading singular DOS correction due to dynami- 
cally screened Coulomb interactionB is shown diagram- 
matically in Fig. ^. It is seen that it is a linear functional 
of the quasiparticle distribution function f{E). It may, 
therefore, be evaluated as a linear superposition of two 
equilibrium electron liquids with respective chemical po- 
tentials /ii — 0, /ifl — ell, and by going over to the 
Matsubara representation. The trace involved in 6N{E) 
may be conveniently evaluated in a basis of momentum 
states p. It reads trSQ = tT{T^Sg^ + r'^6g^)/2T, with 
(Fig. |), 



p 9;f^>l-E| 



vq{in, iE) G'^{iEfG'^_q{iE - in), 

(7) 



where Gp{iuj) = l/{iE + fia — Ep + i/Sr sgnE), a = L,R 
are the non-interacting Green functions. The dynami- 
cally screened CoyJomb interaction in 2D, including ver- 
tex corrections, isEJ 



Vq{Z, Z) 



^v,Q{^z"{z + Z)") 



[Z sgnZ" + iq^D][Z sgnZ'' 



(8) 



with " denoting the imaginary part. In Eq. (g) there 
appears the Fourier transform of the bare Coulomb 
interaction in samples with finite thickness d, Vq = 
{Ae^d/q) arctan(7r/gd), which shows a crossover from 3D 
behavior at large wave numbers (g > 7r/d) to 2D behavior 
at small wave numbers {q < ir/d). The inverse screening 
length in 2D is k = iiie^Nod. According to our sample 
dimensions. In jkp £ ^ d L, the quasiparticle mo- 
mentum p is integrated over the 3D momentum space, 
while the wave number of the diffusion mode is 2D and 
restricted to the range 27r/L < \q\ < 2Tr/£. Thus we 
obtain the DOS correction 



dJ\f{E) = 



1 



EtIi 



1/r 



1/t 



dn , , 1 r 



(9) 



1 + {a^^EthY 

2 + {At:^D/PY 



In 



Vi^ + {2TTKD/Lf 



where Exh = hD/L^ « 0.35 meV is the Thouless energy. 
The first term in quare brackets originates from the q = 
diffusion mode, which must be counted separately due to 
the discreteness of the allowed q values in the 2D sample 
of finite lateral size L. In the case of an infinite, 2D film 
{L oo, Erh < A:BT),|4his expression recovers the well- 
known DOS correctiontl (57V cx -ln(£'r/n) ln{E/hK^D). 



For our finite-size bridges, however, the log divergence of 
the term in square brackets is cut off, and a single loga- 
rithmic singularity remains. Thus, we expect a crossover 
from log^ to simple log behavior in finite-size 2D films, as 
the energy considered becomes less than Eth- Near the 

Fermi edges {\E\ - {2T:fETh, \E + eU\ - {2TrY Eth) the 
logarithmic DOS correction may be cast into the form 



HTt) + 



(10) 



where /(e) = f{Q./kBT) is the distribution function 
in terms of the dimensionless energy. SAf{y, u, T) de- 
pends on two dimensionless energies, the quasiparticle 
energy y = E/ksT and the bias voltage u — eU/kBT. 
At finite T it exhibits a logarithmic dip at each of the 
two Fermi edges at y = and at y = — eU/kBT, as 
shown in Fig. n. The differential conductance correction 
SG/A = [G{U, T) - G(0, T)]/A, with G{U T) = dl/dU = 
d{I/kBT)/du, is calculated using Eqs. (||), ( p^ and (||), 
where the dimensionless quasiparticle energy E/kBT is 
integrated over, yielding a single-parameter scaling form 
in terms of u = U/kBT. The resulting scaling curve does 
not contain any adjustable parameters and agrees quan- 
titatively with the experimental data, as shown in Fig. IJ. 
It is characteristic for log behavior that in Eq. ( p^ ) , and 
consequently also in the conductance G{U,T), the pref- 
actor of the term depending on the dimensionless voltage 
u is independent of T and is equal to the amplitude of the 
T-dependent term In(TT) at y = 0. This is in contrast 
to, e.g., power-law scahng or log^{E/kBT) dependence 
of the DOS correction and, hence, can be used to distin- 
guish the latter behaviors from simple log scaling in the 
experimental data. 

From the above discussion it is seen that the scaling of 
the non-linear conductance in terms of eU/ kBT arises be- 
cause (1) the low-energy behavior of the DOS inside the 
bridge is dominated by one single infrared divergent pro- 
cess (in the present case diffusion modes coupling to the 
single-particle states via dynamically screened Coulomb 
interaction), which ensures that the DOS is described by 
a scale-independent function of the quasiparticle energy; 

(2) the low-energy cutoff, which determines the units of 
the dimensionless energies y and u, is the same in the 
reservoirs and in the interacting region. This means that 
the reservoir temperature T (entering through the dis- 
tribution functions f° in Eq. (^) and the width of the 
Fermi steps in the bridge (entering through Eq. ^) must 
be the same, i.e. energy relaxation in the bridge must be 
negligible; (3) a finite bias voltage is maintained between 
the leads. A high lead-to-bridge aspect ratio and a short 
nanobridge are essential to establish conditions (2) and 

(3) at the same time. 

Finally, we can estimate the amplitude A of the 
Aronov-Al'tshuler conductance correction compared to 
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the background conductance G(0, 1 K). Assuming that 
each channel contributes equally to the conductance, 
the unknown coupling constant T cancels in the ratio 
r = A/G{0, 1 K). We obtain r « 1.0 • 10"'', which is in 
reasonably good agreement with the experimental result 
r ~ 2.6 • 10"'* (see Fig. considering the crudeness of 
the estimate. 




E/eU 

FIG. 6. The normalized DOS correction SAf{E), Eq. (10), 
is shown for symmetric couplings F^ — F^ at fixed tempera- 
ture T for various values of the dimensionless voltage eU /ksT. 

Very recently, it has been proposed that logarithmic 
corrections to the conductance may also be caused by 
charge fluctuations between bridge and the reservoirs 
even in a clean systemEj. This correction can presumably 
be understood as the leading logarithmic behavior of an 
orthogonality catastrophe caused by the vanishing over- 
lap between different charge states of the bridge. A siiH-, 
ilar effect has been considered by Matveev and Larkinc2l 
in the context of tunneling through a charged impurity. 
It remains to be seen whether this may provide an expla- 
nation of the ZBA in our samples. 



IV. CONCLUSION 

We have performed differential conductance measure- 
ments on small 2-dimensional metal bridges, whose 
Thouless energy is larger than the temperature. A con- 
trolled non-equilibrium situation was established by us- 
ing a high aspect ratio between the lead and bridge 
thickness. We observe a ZBA with logarithmic scaling 
behavior as a function of the dimensionless bias volt- 
age eCZ/fesT. The observed scaling curve as well as the 
amplitude of the effect are in good quantitative agree- 
ment with our theoretical prediction that in short bridges 
{Eth > ksT) the Al'tshuler-Aronov DOS correction 



shows Inw behavior instead of the (Inw)^ behavior of 
infinite films. For the observed scaling behavior to oc- 
cur it is crucial that, despite the metallic contact to the 
leads, eigenstates in the bridge survive due to strong 
wave-function mismatch caused by the enormous lead- 
to-bridge aspect ratio. 
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